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Matrix models have wide applications in nuclear theory, condensed matter theory 
and quantum field theory. I discuss supersymmetric extensions of matrix models 
and their applications to branched polymers, the meander problem, and super- 
strings in lower dimensions. 



1 Introduction 

Matrix models have wide applications in nuclear theory, condensed matter 
theory, .high energy theory and quantum gravity since the seminal paper by 
WignerB. An extension to supermatrices can be found in Refs.cHj'Q. nn 
I report in this talk some results on another supersymmetric extension uu 
of matrix models, which is based on one complex bosonic matrix B and fine 
fermionic matrix F and is along the line proposed by Marinari and ParisiQ. I 
discuss applications of the supersymmetric matrix models to branched poly- 
mers, the meander problem, and superstrings in lower dimensions. 

2 Supersymmetric matrix models 

The supersymmetric matrix models in the D — dimensional target space are 
built out of the "superfields" □ 

W a = (B,F), W a = (B\F), (1) 

where a = 1,2 while B and F are general complex bosonic and fermionic 
(i.e. Grassmann valued) N x N matrices, respectively. In other words, the 
hermitean conjugated B^ ^ B and the Grassmann involuted F ^ F. 
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There is no need of introducing a superspace coordinate 9 which would be 
dimcnsionless in D = 0, since the propagators for both bosonic and fermionic 
matrices coincide: 

Bi i B kl) Gausa = Jj S U$kj , ( F H F kl) Gaum = jjSilSkj- ( 2 ) 

Hence, the supersymmetry reduces in D = simply to rotations between the 
B- and i^-components. The proper transformation reads 

8 e B*=Fe, 6 e F = -eB, (3) 

5 s B = eF, S g F = -B j 'e, (4) 

where e and e are Grassmann valued. Note that it is a huge symmetry since 
the parameters e and e are N x N matrices. 

The simplest Gaussian supersymmetric potential reads 

V Gauss = NtvWW, (5) 

where 

2 

WW = ^2 W a W a = B^B + FF, (6) 

a=l 

which reproduces the propagators (Q). It is obviously invariant under the 
rotation (||)-(|]). It is also clear from Eq. (|5|) why one needs complex matrices 
in D = 0: the trace of the square of a fermionic matrix vanishes. 

Any potential, which is symmetrically constructed from the "superfields" 
(Q) , is supersymmetric so that contributions from the loops of the bosonic and 
fermionic matrix fields are mutually cancelled which is the key property of the 
supersymmetry. 

A general interaction potential, which is invariant under the matrix super- 
symmetry transformation (|^)-(^), reads 

^gcn (ww) = NY,jr tr ( Ww ) k ' ( ? ) 

k>l 

where gu are the coupling constants. This invariance can be seen from 

5 e WW = 6 e (B^B + FF) = FeB - FeB = 0. (8) 

The supersymmetric matrix model with the potential (0) describes branched 
polymers as is shown in the next Section. 
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The supersymmetry transformations (Q) can be formalized by intro- 
ducing the (matrix) supercharges ^ 



so that 



(9) 



= [trQe, ...] , S s ... = [ixeQ, ...}. (10) 
Their commutators read 

\Qij 1 Qmn\ {Qij 7 Qmn} 0; (H) 

N ( d d x 



N / d d 

^\ Fik dF^~dK Pkn ) 5mj - (12) 



fe=i 



3 Application to branched polymers 

The partition function of the supersymmetric matrix model with the poten- 
tial (0) reads 

r 2 

Z[g] = / [] dW a dW b e- v ^( ww > =1. (13) 
J 0,6=1 

It is equal to 1 because of the cancellation of bosonic and fcrmionic loops due 
to the supersymmetry. 

Likewise, all the supersymmetric correlators vanish, for example 

/_Ltr (WW) n \ = J dWdWe- y ^( Ww ) ±tr (WW) n = 0. (14) 

Physical quantities of the model are described by the correlators of the 
pure bosonic matrices, which are nontrivial. Their generating function is 



/I 1 \ 1 00 1 

v ' n— 1 



(15) 



with 

G n = (±tr (BlBr}. (16) 

''Here and below the order reflects matrix multiplication. 
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Figure 1: Typical diagrams for Gi when the potential yjlfW) is cubic in WW (above) 
and the associated branched polymer graphs (below). 



The imaginary part of G(A) determines the distribution of eigenvalues of the 
matrix B and, therefore, the spectrum of the proper statistical model. 

The correlators (16) can be calculated □ using the Schwinger-Dyson equa- 



tions and the supersymmetry Ward identities. As a result, G(A) obeys at large 
N the closed equation 

(AG(A) - 1) V (A - 1/G(A)) = AG 2 (A) , (17) 

whose limit of A — * oo yields the equation 

G x V'{Gx) = 1 (18) 

for G\. This equation is quadratic for a quartic potential. 

The fact that a closed equation is obtained for the propagator Gi is a 
consequence of the cancellations between bosonic and fermionic loops. Some 
of the diagrams which survive the cancellation are shown in Fig. ^. For obvious 
reasons, they are called the "cactus diagrams". Note that the diagrams have 
orientation: the cactus loops can only proliferate on the exterior of already 
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existing loops. This is in contradistinction to related bubble diagrams one 
encounters in the large- N limit of pure bosonic or fermionic vector models. 

The critical behavior of the model arises when Eq. (|l^) holds simultane- 
ously with 

[G C V'(G C )}' = ... = [G c V'(G c )} (m - 1] = 0, (19) 

which can always be achieved by tuning m couplings gk of the potential ([?]) . 
Near the critical point, G\ behaves as 

d ~ G c - const • (a c - a) 1/m , (20) 

where a stands for an overall scale of g^s. The susceptibility \ = dGi/da at 
the critical point scales as 

X ~ (ac - «)~ 7str (21) 

with 7 s t r = 1 — which coincides with the (multi-)critical index of the 

branched polymers u. 

This relation to branched polymers becomes explicit by noting that the 
cactus graphs which survive the supersymmetric cancellation have an inter- 
pretation as branched polymers, with the couplings —g^ associated with the 
branching weights. Cutting each loop in a succession, such that only the stat- 
ing line is attached to the vertex, produces a branched polymer graph of a 
"chiral" type since branching only occurs at one side of the open line, corre- 
sponding to the fact that the cactus loops can only be attached to the exterior 
of already existing loops. This is illustrated in Fig. [I]. Equation ( |l8| ) can then 
be rederived pure combinatorially. 



4 Application to the meander problem 

The meander problem is to calculate combinatorial numbers associated with 
the crossings of an infinite river (Meander) and a closed road by 2n bridges]^ 
Neither the river nor the road intersects with itself. These meander numbers, 
M n , obviously describe the number of different foldings of a closed strip of 2n 
stamps or of a closed polymer chain. 

The generating function of the meander numbers can be. represented via 
the following correlator in the supersymmetric matrix model 



M(c) ee V c 2n M n = lim /— tr BB^ In ( f dfadfo e " 

^-J JV->oo \N \J 



, (22) 

Gauss 



c See Ref.LJ for an introduction to the subject. 
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where (f>\ and </>2 are N x N hermitean matrices, the action S is given by 



S = y tr $ + — tr <$>\ - cN tr (faB^ <j> 2 B) - cN tr (faFfoF) (23) 

and the Gaussian averaging is with respect to the action (|J). The presence 
of the log in Eq. ( p2| ) leaves only one loop of the field 4> associated with the 
road while the supersymmetry kills the loops of the field W associated with 
the river. The limit of large N is needed to keep only the planar graphs as in 
the original meander problem. 

Expanding in c, the coupling constant of the quartic interaction, the me- 
ander numbers can be represented as the sum over words built out of two 
letters: 

2 / i 

Mn = V ( 1 r T ^BW a2 ---W a2n _ 1 W ( 



Gauss 



^tvW a2n W a2n _ 1 ---W a2 B) , (24) 

iV / Gauss 

where the order of matrices is essential for the fermionic components. Equa- 
tion ( p4| ) is a nice representation of the meander numbers which looks more 
natural than the one based on the replica trick in a pure bosonic model. 

Equation ( p4[ ) can be represented in an alternative form by introducing 
noncommutative variables u,v and u\v^ which are annihilation and creation 
operators in a Hilbert space with the vacuum |fi) and obey the Cuntz algebra 

uu' = 1, vv^ — 1, uv^ — 0, vv) — 0, (25) 

as well as the completeness condition 

u^u + vh = i-\n) (n\ . (26) 

There are no more relations between the noncommutative variables. 
Denoting 

u a = (u,v), u a = (u,-v), (27) 

the generating function (^2| ) can alternatively be represented as the vacuum 
expectation value 

M( C ) = (niGu^uG^n) = - (n|GwW|n) , (28) 
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where G is given by the continued fraction 
G{u) = L_ . (29) 

1 - y/cu ai j U ai 

1 \/~CU a2 J ^tt2 

1 - \/cu a3 j u„3 

1 - VCM Q4 - U ai 



Here u and u interchange in the consequent lines. The following notations are 
used in Eq. (§§): 

G = G(u), G = G(u), G^ = G (V) , G t = G(M t ). (30) 



The two expressions on the right hand side of Eq. ([28|) are equal due to the 
supersymmetry. 



Though Eq. (28) reproduces the meander numbers when expanded in c, it 
seems to look like a reformulation rather than a solution to the problem since 
it is not clear how to deal with functions of noncommutative variables. 



5 Application to superstrings? 

The critical index of the string susceptibility 7 str for a superstring embedded 
in a ZJ-dimensional space had been calculated from the super-Liouville theory 
and reads Ii2l 

D-l-^/(l-D)(9-D) 

7str = 1 • (31) 

Many (not yet successful) attempts of discretizing. superstring are per- 
formed starting fromllil. A progress has been achievedl!3 only for the simplest 
case of pure 2-dimensional supergravity which can be associated with a su- 
pereigenvalue model. It reveals the super- Virasoro algebra associated with the 
Neveu-Schwarz sector of the superstring. 

The idea to verify whether or not a super- Virasoro algebra can be realized 
in the supermatrix models is to construct the matrix generators (cf. (Q), (p2[)) 

N ( d d \ N ( d d 

Here C%j is Grassmann even and Gij is Grassmann odd. 



(32) 
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The operators Cij and Gij obey the commutation relations 

l/^ij ) ^mn] ^inCmj Ci n S rn j ) (33) 

[^27 : -^mn] ^inGmj Gin^mj : (34) 

\Gij 1 Gran} — ^inC'/nj ~l~ Ci n S m j • (*^) 

These can be derived by explicitly commuting the operators (|3 
The commutator (B3h itself implies the Virasoro algebra 



[C s , Ct] = (t- s)C s+t , C s = tr (£ (WW)") (36) 

for s,i>0asiV-^oo. Given (|36|), the potential (0) can then be recovered 
from the Virasoro constraints 

= J dWdWC s e ( Ww ) = L s [g] Z [g] . (37) 

The conjecture is that the whole matrix algebra (|33])~(|3"5|) implies, as 
JV — ► oo, the super- Virasoro algebra associated with the Ramond sector of 
the superstring in the D — dimensional target space: 

[C s , C t ] = (t-s)C s+t , [G r , C s ] = (| - rj Gr+s , {G r , G s } = 2C r+s . (38) 
An explicit form of the operators C s and Gr can be constructed starting from 
£ =tr£, g =trG, C x = atxCWW + btrQWAW (39) 

where 

2 

WNW = e ab W a W b = B^F - FB (40) 

a,b=l 

and using d38|). For example, one gets 

Gi = ctrGWW + dtrCWAW, 

C 2 = c(c + d)trC(WW) 2 + c(c~d)trGWAWWW (41) 

with 

c = -26, d = 2a + 46 (42) 

and so on. The constant b ^ in C\ since Ci vanishes otherwise. In order 
for this procedure to be successful, all the operators C s and Gr are to be 
nonvanishing. 

The action of the proper supersymmetric matrix model should then be 
determined to reproduce these super- Virasoro operators. It should involve 
both WW and WNW given respectively by Eqs. (|^) and ( f40| ) and, therefore, 
Grassmann odd coupling constants in addition to those in (|7|). 
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